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SUMMARY 

When a closed body or a duct envelope moves through the atmosphere, 
air pressure and temperature rises occur ahead of the body or, under ram 
conditions, within the duct. If cloud water droplets are encountered, 
droplet evaporation will result because of the air-temperature rise and 
the relative velocity between the droplet and stagnating air. It is 
shown that the solution of the steady-state psychrometric equation pro- 
vides evaporation rates which are the maximum possible when droplets are 
entrained in air moving along stagnation lines under such conditions. 
Calculations are made for a wide variety of water droplet di ame ters, 
ambient conditions, and flight Mach numbers. Droplet diameter, body 
size, and Mach number effects are found to predominate, whereas wide 
variation in ambient conditions are of relatively small significance in 
the determination of evaporation rates. 

The results are essentially exact for the case of movement of drop- 
lets having diameters smaller than about 30 microns along relatively long 
ducts (length at least several feet) or toward large obstacles (wings), 
since disequilibrium effects are then of little significance. Mass losses 
in the case of movement within ducts will often be significant fractions 
(one-fifth to one-half) of original droplet masses, while very small drop- 
lets within ducts will often disappear even though the entraining air is 
not fully stagnated. Wing-approach evaporation losses will usually be of 
the order of several percent of original droplet masses. 

Two numerical examples are given of the determination of local evap- 
oration rates and total mass losses in cases involving cloud droplets 
approaching circular cylinders along stagnation lines. The cylinders 
chosen were of 3.95-inch (10.0+ cm) diameter and 39.5-inch (100+ cm) 
diameter. The smaller is representative of icing-rate measurement cylin- 
ders, while with the larger will be associated an air -flow field similar 
to that ahead of an airfoil having a leading-edge radius comparable with 
that of the cylinder. It is found that the losses are less than 5 per- 
cent. It is concluded that such losses are, in general, very small 
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(less than 1 percent) in the case of smaller obstacles (of icing-rate- 
measurement -cylinder size); the motional dynamics are such, however, that 
exceptions will occur by reason of failure of very small droplets (moving 
along stagnation lines) to impinge upon obstacle surfaces. In such 
cases, the droplets will evaporate completely. 


INTRODUCTION 

Several studies (for example, refs. 1, 2, and 3) have recently been 
made of trajectories of water droplets entrained in atmospheric air 
approaching wings or liquid -water-content measurement cylinders. In 
these trajectory studies it has been assumed that the droplet shape and 
mass remain unchanged despite possibilities of deformation and evapora- 
tion arising in connection with the presence of local compressible flow 
fields in the vicinity of wing or cylinder. In such fields, air temper- 
atures and pressures are no longer those of the undisturbed stream. 
Moreover, the droplet acquires, in general, a velocity with respect to 
the disturbed air. 

In view of the air temperature and pressure changes that occur and 
the existence of an air-droplet relative velocity, it appeared desirable 
to determine the magnitudes of evaporation mass losses occurri ng in the ' 
vicinity of obstacles or within ducts. The assumption that the droplet 
remains spherical was retained. 

The general literature in the field of evaporation from droplets is 
extensive, but few studies of the particular problem of droplets approach- 
ing obstacle under icing conditions have been made. Both Hardy (ref. 4) 
and Langmuir (ref. 5) have considered certain aspects of the question; 
both analyses took into account the fact that the droplet will not, in 
general, be in instantaneous psychrometric equilibrium with its surround- 
ings . In both cases, the conclusion was reached that total evaporative 
losses from droplets actually reaching obstacle surfaces are of the order 
of several percent for the particular sets of conditions of their analy- 
ses. In both investigations, however, the ranges of droplet size, air 
temperature and pressure, body size, and flight Mach number were rather 
restricted. 

A complete treatment of the problem would require: 

(1) use of the equations of motion of a droplet of liquid in a (com- 
pressible) gas, the equations to take into account in some way changes of 
drag coefficient with Mach number; 

(2) use of the equations governing the dynamic thermal behavior of 
a volatile sphere under quite general conditions of changing heat and 
mass-transfer rate including, in some instances, radiation effects; and 
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(3) coverage of wide ranges of obstacle size and shape, droplet 
size, flight Mach number, and flight ambient air conditions. 

The unavailability of a suitable high-speed calculator made it 
necessary, however, to restrict in some manner the scope of the calcula- 
tions. 

Since the subject of droplet evaporation, insofar as aircraft are 
concerned, is of interest chiefly in connection with icing, radiation 
effects could be ignored (since all temperature levels involved are low). 

Further, it was decided to consider that a droplet remains in in- 
stantaneous psychrometric equilibrium with its immediate surroundings, 
and its internal temperature was to be taken as uniform. Certain quanti- 
tative aspects of this assumption are scrutinized in the body of the re- 
port and in appendix B; a qualitative justification of the procedure is, 
however > given here. 

If evaporation occurs at all, it will occur principally as a result 
of a droplet temperature rise, with which a droplet -surface vapor -pres sure 
rise will be concomitant. Therefore, maximum evaporation rates will occur, 
in virtually all cases, in the vicinities of body stagnation points. If, 
now, bodies and flight conditions (wing attitudes, for example) are con- 
sidered such that air flows in the vicinities of stagnation points are 
essentially symmetrical about such points, then the trajectory of a drop- 
let originally on a stagnation line will essentially coincide with the . 
line. For the sake of simplicity, then, let attention be confined to 
histories of droplets moving along such lines. Air temperature and pres- 
sure will rise monotonically. The droplet instantaneous position will 
always be ahead of that of the air with which, at an arbitrary previous 
time, it was in contact, although a quasi -static mass-transfer analysis 
does not demand such motion. More importantly, with or without evapora- 
tion the interior droplet temperature will always lag behind the rising 
surface temperature, and both the mean droplet temperature and the sur- 
face temperature will always be less than the local equilibrium (psychro- 
metric) value; all droplet temperatures will rise monotonically. It is 
possible to conclude that both the acual local time rate of evaporation 
and the actual total loss of liquid will be less than those calculated 
on the basis of psychrometric calculations in the case of stagnation 
streamline droplet motion. It therefore follows that the present quasi - 
static calculations set upper bounds to the loss rates and total losses 
for motion along stagnation lines. 

Finally, evaporative losses off stagnation streamlines will in gen- 
eral be less than stagnation-streamline values because of the fact that 
maximum temperature rises occur along the latter. 

Accordingly, it is proper to consider that psychrometric-equilibrium 
calculations establish limiting local (and, when integrated, total) losses 
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which are never exceeded in real situations. Moreover, the results 
apply without change to the histories of small droplets moving within 
long ducts or toward large obstacles. 

Since the calculations of mass losses are made under the assumption 
of the existence of successive instantaneous psychrometric -equilibrium 
states, it is desirable that the degree of quantitative applicability of 
that procedure be ascertainable in any given case. In the section of 
the report entitled Response of Sphere to Changes of Environmental Tem- 
perature, two parameters are mentioned by the use of which it is possible 
to assess in a rough manner the degree to which such equilibrium assump- 
tion calculations lead to values concordant with mass-loss rates which 
would be obtained in the case of more rigorous calculations. The expres- 
sions determining the magnitudes of the two parameters are derived in 
appendix B. All symbols used are defined in appendix A. 

The detailed calculation procedure used to obtain the mass-loss 
rates themselves is presented in the section of the report entitled 
EVAPORATION CALCULATIONS. 

The results and possible applications are then discussed; two numer- 
ical examples are given. All symbols used are defined in appendix A. 


General Considerations 

The assumption that instantaneous local psychrometric equilibrium 
exists between droplet and atmosphere eliminates dependency of calcu- 
lated local loss rate upon previous droplet history. This is true be- 
cause local droplet calculated temperatures then depend only upon in- 
stantaneous local conditions. 

When that assumption is made, calculations of maximum possible total 
mass losses in actual physical situations may be made in the following 
manner : 

Droplet mass-loss rates for reasonable combinations of droplet size, 
local ambient (relative) conditions, and relative air -droplet speed are 
calculated; the details of the procedure are given in the section en- 
titled EVAPORATION CALCULATIONS. 

The variations of air speed, temperature, and pressure along the 
forward stagnation line are obtained in some manner for the particular 
physical situation of interest. 

The motional history of a droplet of the size for which the calcula- 
tions are being made is then obtained for the case of motion essentially 
along the stagnation line; this was accomplished by the use of a differ- 
ential analyzer or by the use of some numerical technique. The local 
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rates of evaporation were then obtained simply by interpolating among 
the results of the generalized evaporation calculations at reasonable 
intervals; the tacit additional assumption is made here that the over- 
all mass loss will be small, so that constancy of droplet size may be 
assumed to a sufficiently high degree of approximation in making a cal- 
culation of local mass -loss rate. 

Finally, the over-all mass loss is obtained by adding the small loss 
increments corresponding to the selected droplet path increments. If 
the over -all loss is greater than, say, 10 percent, the decrease of drop- 
let size can easily be taken into account. 


RESPONSE OF SPHERE TO CHANGES OF ENVIRONMENTAL TEMPERATURE 

It has been found practicable, in the calculations made in connec- 
tion with the present study, to cover wide ranges of the several variables 
involved only because droplet internal temperature uniformity and psychro- 
metric equilibrium with the environment were assumed. These assumptions 
greatly simplified the calculations while at the same time permitting 
conclusions to be reached as to limiting evaporation rates in real situa- 
tions . 

It was considered desirable, however, to attempt to establish cri- 
teria whereby the degree of departure of equilibrium values of droplet 
temperature and evaporation rate from actual values could be estimated. 
While precisely this goal was not attained, criteria were established 
which serve to determine what might be designated the "degree of invalid- 
ity" of the basic equilibrium assumptions. 

The criteria themselves are derived in outline form and are discussed 
in appendix B. It is sufficient here to mention that three variables are 
of importance in this connection. Two of these are the times required for 
a sphere to respond to sudden surface and to environmental changes of 
temperature. The third variable is the "approach period", which is arbi- 
trarily taken as the period required for the droplet, moving at a speed 
equal to the relative air -obstacle free-stream speed, to move through an 
approach distance equal to a principal transverse dimension of the ob- 
stacle. At the latter distance upstream of the stagnation point, the 
air -flow characteristics, while still substantially those of the free 
stream, will have departed by easily measurable fractions from those of 
the free stream. For example, the stream-obstacle relative Mach number, 
for a free-stream Mach number of 0.75, is 0.74 at a distance of one 
radius from a cylinder surface. It follows that the approach period, as 
here defined, represents in approximate fashion the time during which 
most of the change of effective droplet environment occurs as the drop- 
let approaches the obstacle. 
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The two parameters which then determine the degree of invalidity of 
a quasi-static treatment are any two among the following three: Ratio 

of internal to external sphere time constant (for response to a sudden 
surface or environment change, respectively), ratio of internal time con- 
stant to approach period, and ratio of external time constant to approach 
period. 

In general, if both the ratios of time constant to approach period 
are small, a quasi-static (psychrometric-equilibrium, uniform internal 
sphere temperature) analysis is quantitatively exact. In all other cases, 
as indicated in appendix B, the situation is complex and does not read- 
ily yield to any simple analysis. 

In particular) a quasi-static analysis is exact in the case of the 
approach of smaller droplets (below about 30 microns diameter) to large 
obstacles (wings) or motion of such droplets within ducts at least sev- 
eral feet in length. In the case of the long duct, in which, except for 
a region near the mouth, conditions are uniform, total evaporation losses 
are obtainable directly from the results of calculations made in connec- 
tion with the present study, as indicated in the section entitled RESULTS 
AND DISCUSSION. 


EVAPORATION CALCULATIONS 

The basic numerical result required in these studies was the rate 
of loss of mass of a water droplet of a particular size moving at a pre- 
scribed Mach number with respect to air characterized by a specified de- 
celeration from the original free -stream flight speed and by the free- 
stream temperature and pressure. 

As discussed previously, it was assumed that the droplet remained 
in instantaneous psychrometric equilibrium with its effective environ- 
ment; the effective environment consisted of the mixture of air and water 
vapor the local apparent temperature and pressure of which, with respect 
to the droplet, were assumed to be altered from the true local (static) 
values by the motion of the droplet with respect to the air. The quanti- 
tative nature of the assumed changes is discussed subsequently. 

In outline form, the calculation scheme was the following: 

Ambient static air temperature T gt *, ambient static air pressure 
(inclusive of water vapor P gt ^J, and flight Mach number were 

assumed. A deceleration of the air (fig. l) to a selected relative air- 
obstacle Mach number Mr was then assumed, and the increased local air- 

plus -water -vapor static tender atures and pressures were computed; the 
change was assumed adiabatic and isentropic. The effective air and 
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water-vapor temperature and pressures, on the basis of a given air- 
droplet relative Mach number M r , were then calculated. 

It was then possible to compute a local relative air-droplet Reynolds 
number Re r , from which heat- and mass-transfer coefficients were ob- 
tained. 


By using two temperature -range overlapping quadratic approximations 
for the vapor pressure of water as substitutes for the more accurate ex- 
ponential expression, it became possible at this point to obtain, begin- 
ning with the psychrometric (he at -balance) equation, a quadratic equation 
in the droplet temperature and other quantities. This equation was then 
solved for the droplet temperature; the procedure avoids the iterative 
calculations otherwise required. The computed mass-transfer coefficient 
and droplet vapor pressure were then used to obtain the required rates 
of evaporation. 

The significant steps of the calculation are given herei naf ter with 
appropriate explanations of symbol meanings and origins of relations 
where those are not immediately obvious; the fundamental relation is the 
psychrometric equation 

^d,e,Z " ^d,s^h = ^ tl m^Pd,s " p d,e,Z^ (l) 

which indicates that the droplet temperature assumed a value such that 
as much heat is removed in unit time by evaporation as is supplied by 
free or forced convective heat-transfer processes. 

The initial steps are required to determine the effective tempera- 
ture, air pressure, and vapor pressure of the local environment with re- 
spect to the droplet. 

The static temperature at the instantaneous droplet position (here- 
after identified as local) is a function of the air-obstacle relative 
total temperature and the local air-obstacle relative Mach number: 



1 + 0.2 ^ 
1 + 0.2 Mg 


The mean effective temperature of the environment with respect to 
the entire droplet is then calculated from the local static temperature 
and the droplet -air relative Mach number: 


T d,e,Z = T st,Z d + 0-16 M^) 
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Note that a recovery coefficient of 0.8 is used; this is no more than 
a representative average figure based on many data reported in the liter- 
ature. 

Implicit in this relation is the concept that the laws governing 
events occurring at the forward position of the sphere may be applied 
with little error to the entire sphere. Fortunately, the Reynolds num- 
ber regime (0— * 1000) of interest is precisely that in which, for most 
smooth, regularly shaped bodies, heat losses and mass -transfer rates are 
much greater over the forward portions of the bodies than over the rear 
(unpublished interferometric research performed by E. R. G. Eckert at 
Wright Field indicating that, in the Reynolds number range 24.5$Re^605, 
front-to-back local Nusselt number ratios are of the order of 3 or 5 
to 1; see also refs. 6 and 7). It is therefore permissible to apply the 
law of variation (of the particular quantity under consideration) for the 
forward portion to the entire body, although the absolute level (for ex- 
ample, of a heat-transfer coefficient) will be determined by an experiment 
involving the whole body at some typical Reynolds number . 

Insofar as temperature effects are concerned, there will of course 
be some heat conduction from the front to the rear, but the fact that the 
film resistance is much higher in the rear will ensure attainment of 
essentially front-surface effective temperatures . 

The local air-droplet relative total pressure is used in the calcu- 
lation of air-droplet relative Reynolds numbers, and the relative total 
vapor pressure in the calculation of air-droplet relative vapor pres- 
sure. The reason for this, apart from the fact that little evaporation 
occurs over the rear portion of the droplet, is that it is shown in ref- 
erences 8 (fig. 8) and 9 (table IV, fig. 3) that at Reynolds numbers be- 
low about 50 the pressure in the vicinity of the forward stagnation point 
is substantially higher than " total" . 

Emphasis is placed on the low Reynolds numbers regime because air- 
droplet relative Reynolds numbers will, generally, be small for the 
smaller (say, below 20 microns) droplets, while evaporation rates will 
be highest for such droplets. 

Accordingly, a rough average pressure for the front surface is the 
total pressure; the final results, in any case, are not critically de- 
pendent upon the nature of the laws of pressure variation assumed in 
these calculations. 

The local air-droplet relative total pressure is given by 

p ( 1 + 0.2 M 2 ) 7 / 2 (1 + 0.2 M 2 ) 7 / 2 

e,7,r st,oo 


(1 + 0.2 M 2 ) 


2 ^ 7/2 
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In the calculation of effective local water vapor pressure, it is 
first necessary to determine the original (static, free-stream) vapor 
pressure , either complete or 90-percent saturation being assumed. 

While it is necessary only to consult a table of values to deter- 
mine the vapor pressure in question, it was decided that the quadratic 
relations used subsequently t'o calculate the saturation vapor pressure 
at the droplet surface should also be used in other vapor-pressure cal- 
culations so as to avoid the small inconsistencies which would otherwise 
arise. 


The vapor pressure in the free stream is given by 


p, = c( 6136 + 438.7 t . + 10.69 t . J 

d ,oo st,eo st , 0 # 


for 


-15° < t < 15° C 

st,» 


and 


p, = a (5288 + 287.7 t . + 4.3047 

d ,oo st, « 


t 2 + ) 

st /» 


for 


(2a) 


(2b) 


-30° C < t st)B < -8° C 

These relations yield approximately correct values of vapor pres- 
sure over the intervals indicated. The accuracy of the second equation 
is higher over the common interval -15° C-* -8° C; the errors at var- 
ious temperatures are listed in table I (for the case a = l). The 
accepted figures were taken from reference 10. 

Although neither of the two vapor pressure - temperature relations 
used is accurate above 15° C, the first one given was used for droplet 
temperature calculations for droplet temperatures as high as 25° C; the 
results are merely qualitative in the range 15° to 25° C, but are re- 
tained for the sake of completeness. 

In the calculation of effective local ambient vapor pressure, the 
assumptions are made that the ratio of vapor to air pressure remains 
constant and that y for a mixture of water vapor and air is essenti all y 

1.4. 


In connection with the calculation of droplet-air relative Reynolds 
number, a temperature at which the viscosity is to be evaluated must be 
selected. The droplet-air relative total temperature 
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(1 + 0.2 M ^) 

T = T — 

d,t,l d,e,l + 0>16 M 2) 

was selected for that purpose in the earlier stages of the work and sub- 
sequently retained for the sake of consistency despite the fact that vise 
of a relative total temperature cannot be justified. The maximum dif- 
ference between the relative total and relative effective temperatures, 
however, occurs at Mj. = 0.75 (an unusually high figure) and is ~5° C. 

The difference between the viscosities corresponding to the two temper- 
atures is then "1.5 percent, while for Mj. values below 0.5 the dif- 
ference is negligible. Since the viscosity is used only to calculate a 
droplet-air relative Reynolds number, the effect of any such difference 
is virtually undetectable . 

The droplet-air relative Reynolds number is given by the expression 

p. _ 2p st,Z r d v r 

ne r 


which is, accordingly, a Reynolds number based upon local static air 
density, droplet diameter, relative air -droplet speed, and viscosity at 
^d,t,7* ^ le viscosity, to avoid irregularities resulting from the use 

of tabulated experimental values, was obtained from the expression 

\i (poise) = 10 -6 [l72.81 + 0.48721t - 0.0004932 t 2 + 9.33X10 -7 t 3 ] 

The relation given subsequently between the Nusselt number for heat 
transfer and the Reynolds number is based primarily on the survey made 
by Williams (ref. ll). A Prandtl number of 0.73 was used in the calcu- 
lations. The exponent 0.6393 was selected, in part, because it had pre- 
viously been found that the drag., coefficient of a sphere is given quite 
accurately over the range 0 < Re < 800 by a relation in which that 
exponent appears; it had therefore been possible, by using that exponent 
for the heat-transfer equation, to save effort in certain calculations. 
The expression fits the experimental data well over the range 
0 < Re < 10 4 . 

Nu^ .= 2.000 + 0.2464 Pr 1 / 3 Re 0.6393 

It should be noted that the assumption has been made, in connection 
with the expressions for the heat-transfer Nusselt number just given and 
for the mass-transfer Nusselt number given later, that heat- and mass- 
transfer processes do not affect each other (except for the establish- 
ment of the psychrometric temperature). Otherwise stated, the assumption 
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is made that the presence of water vapor in the sphere boundary layer 
does not affect either the inert gas (air) properties or the eddy diffu- 
sivity and thermal conductivity. While it is believed that this assump- 
tion is valid at low air temperatures and (in particular) at low rates 
of evaporation, evidence exists tending to invalidate the assumption 
under other conditions (ref. 12). 

The Nusselt number for mass transfer subsequently is somewhat less 
familiar than Nu^; it is defined as (hjjjD^/k^ jin ) , where is the 

mass-transfer coefficient in cm"^ sec”-'- (dyne cm'^)' 1 and kj^ is the 
mass-transfer "conductivity" defined as (p^/F a ) . In the latter, 

p is the effective medium density, g the diffusion coefficient of the 
pair of substances, and P a a local mean particle pressure of the gas 

into which the diffusing material is passing (in this case, dry air). 

WUjh = 2.000 + 0.330 ReJ- 56 

The coefficient 0.330 is actually the product of the constant 0.39 
and of the selected value of the one -third power of the Schmidt number 
(0.605) in the expression: 

Nu = 2.000 + 0.39 Sc 1 / 3 Re 0 * 56 
m r 


The exponent of this relation is that of Williams (ref. ll), while the 
constant was selected to ensure good fit over the range 0 <Re < 3000. 

It will be noted that Williams ' data are plotted in terms of the Colburn 
parameter j m , which is defined as follows (notes m = (gp a /P a ) (M /m )) 



= Ru Re' 1 
m 


pg ^a ^d\ 

Pa p 2/3 ^1/3 p 2/3 



= 1% Re' 1 
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It therefore follows that Nusselt numbers of mass transfer are obtain- 
able from Williams ' (and similar) data by multiplying the values 

by Re Sc 1 / 3 . Heat -transfer data published since 1942 agree reasonably 
well with the present empirical equation used herein; mass-transfer data 
exhibit greater scatter. 

It is possible to select a fixed value of Schmidt number because 
the latter varies rather slowly with both temperature and pressure for 
rather wide ranges of temperature and pressure in the vicinity of stand- 
ard atmospheric conditions. It is possible to deduce this slow varia- 
tion ■upon inspection of the relation 

* ■ yy •" 

of page XVI -16 of reference 13. (Pq is 0.220 cm 2 /sec at 0° C for the 
system water -air, ref. 14.) Since, in the expression Sc = h a /p a P> the 

viscosity varies roughly as the 0.8 power of the temperature, while air 
density varies inversely with temperature and directly with pressure, 
the residual variation of Sc is a rather slow decrease with tempera- 
ture. 


The psychrometric equation (l), for convenience in discussing the 
remaining calculations, is how restated in the form 


^s " ^e,! " x (Pd,s ' p d,e,l) 


where 


Ah Ak Nu 

^ _ m __ i ,m m 

" h h *f,h Nu h 

The thermal conductivity k_^ ^ is nearly independent of pressure at 
moderate pressures, whereas kf >m may be written 

*f,m - (N.P./P) 

In this instant, P = P e ^ r » The quantity X may now be expressed as 
follows : 


X = 


_ /V 


nr 


Nu h \ k f, h/ 


N.P. 


1,013,250’ 
P , 

^ e,Z,r / 
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Further, the quantity (kf >m Af ^N.P. increases only about 0.11 percent 

per °C increase of temperature at 0° C. For present purposes, therefore, 
the expression 

X(N.P.) p ^ fixed temperature 

was evaluated once at a temperature in the vicinity of -8° C from the 
following relations: 


X = 2.519X10 10 erg/g 
[kf,jJ-8° C = 2422 erg/cm sec °C 

[j^f ,m] N.P. ,-8° C = 1-746X10- 10 g cm 2 /erg sec 
The result was (543.5X10-0)” 1 , so that, finally. 


1 c Nu h 

r - - 1 - = rar RXin- b p . 

Nu m 


= 543.5Xio _D P„ 7 „ 

e , 0 f 1 


Constants Cj, C 2 , and C 3 are now defined by the relation 

2 


p a, £ 


“ C 1 + C 2 t d,s + C 3 t d,s 


(3) 


where the constants have been given previously in the quadratic expres- 
sions (eq. (2)) for vapor pressure. Equation (3) and the psychrometric 
equation (eq. (l)) are then used to obtain 


t2 + 
d,s 


1 + X C 2 
XC* 


C 1 - Pd,e,Z _ ^d,e , l 


d,s 


XC, 


which may also be written 


In the latter relation 


t , + a., t , + a_ = 0 

d,s 1 d,s 0 


a l - S 1 ^ + °2> 


and 


a n = C"-*- (C, - p, 7 - t , 7 X - - 1 -) 

0 3 1 *a,e,Z d,e,Z 



14 


NACA TN 3024 


(The choice as to the particular set of constants to he used is made on 
the basis of a guess as to the final value of t^_ s itself in border- 
line cases.) 


The droplet temperature is then given by the relation 



The rate of loss of mass is obtained from the expression 


Qm “ mA d = kjjjAd ( p d,s * p d,e,l) 

in which m is the mass-loss rate in g/cm 2 * * * sec, while A d is the 
surface area. Now, 

^f,m ^ u m 

“"■-'ST - 


and 



(kf^N.P.T. 


/T \°' 89 
1,013,250 ( 1 d,e,r) 

p e,Z,r \273.2 / 


From these. 


% - 2jtr d( p d,s-Pd,e,z) 


( k f,m)N.T.P. (l,015,250)Nu m 


e,Z,r 


/ T d,e,r\ / 265.2 V 
\273.2 J \273-2 J 


- 0.11 


is obtained; a factor (T d e r /273.2) ” 9 ’ ^ has here been replaced by 
- 0*11 } y 

(265.2/273.2)” ’ to facilitate calculations. The procedure is 
admissible because of the negligible variation of the factor over the 
temperature range of interest. The numerical constant then equals 


2 jt XI. 746X10 10 (0.9707) ”°* 11 x = 4.082X10 -6 

and the following expression yields the desired mass-loss rate: 

4.082X10” 6 r d T a>ejr (p djS - P d , e ,z) Nu m 

Qm “ P 

e,Z,r 
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The time rate of loss of fractional mass is then given by 


A = 


30m 

4jrpdr^ 


The distance rate of loss of fractional mass may be calculated in an 
elementary manner only if the droplet speed is always positive with 
respect to the air, and the droplet trajectory and an air streamline 
virtually coincide. If these conditions are met, the rate Q is given 

fcy 

Cl = V( u l + V r ) 

In tabulations and graphs Q, which as calculated has the dimensions 

cm -1 , appears as the fractional mass lost per foot of droplet travel 
toward the obstacle. 


RESULTS AND DISCUSSION 

The calculations were made for a number of combinations of ambient 
relative humidity (unity and 0.9), flight altitude (3, 6, 10, 15, 22.5, 
and 30 thousand ft), ambient temperature ( -30, -25, -20, -15, -10, and 
-5° C), droplet diameter (5, 10, 15, 25, 50, and 150 microns), flight 
Mach number (0.25, 0-50, 0.75, and 1.0), air-obstacle relative Mach 
number (0, 0*125, 0.25, 0.50, and 0.75), and droplet -air relative Mach 
number (0, 0.063, 0-125, 0.25, 0.50, and 0.75). Not all of the possible 
combinations were used, however, since the number of such combinations 
is prohibitively large. 

The combinations for which calculations were actually made are 
indicated in the table (table II) in which the results are also pre- 
sented. The air -obstacle and air -droplet results are not tabulated 
because they are contained in their entireties in the flight Mach num- 
ber results collocation (table Il(a) and 11(b)). Table III contains 
a summary of certain entry data (air pressures at various altitudes and 
temperatures) used in the calculations. Most of the results, with the 
chief exception of the entire 90 percent humidity group, are also 
presented in graphic form in figures 2 to 9. 

Figure 2 exhibits the fractional mass lost per foot of droplet 
travel (along a stagnation line) as a function of flight Mach number M» 
at a succession of fixed droplet-air relative Mach numbers M r for sev- 
eral local air-obstacle Mach numbers Mr under ambient conditions as 
indicated (altitude 10,000 ft, air temperature, -25° C), the droplet 
diameter remaining at 15 microns. It will be noted that the evaporation 
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rates increase sharply with increasing M®, particularly at the higher 
air-obstacle relative Mach numbers. For example, at an % of 0 and 
an M r of. 0-125, ft changes from 0.0047 per foot at an M® of 0.25 to 
0-061 per foot (a factor of 13) at an M a of 0.75; at an % of 0.50 
and an % of 0.125, ft changes from 0.00014 per foot at an M» of 0-5 
to 0.0051 (a factor of 36) at an M® of 0.75* 

This relatively large increase of evaporation rates (on a unit 
distance basis) is the result of combined action of at least two effects. 
First, at the lower % values (that is, at nearly full stagnation), 
the time rate of evaporation will tend to be high irrespective of the 
Moo value, although it will, of course, increase fairly rapidly with 
increase of the latter. Second, the increase of droplet speed (with 
respect to the obstacle) at increasing Moo tends to offset the increase 
of time rate of evaporation insofar as the fractional loss per unit 
distance is concerned. The net result is a moderate rate of increase of 
fractional loss per unit distance with increasing Moo at lower Mr 
values. At higher % values (greater than, say, 0-5), however, the 
time rates of evaporation at Mqo values such that the Mr values in 
question involve o nly relatively slight stagnation will be very small. 
Relatively small increases of Mqo, therefore, will markedly increase 
what may be called the fractions of full stagnation, and there will be a 
substantial increase of time rates of evaporation under those circum- 
stances. Because the (subsonic) M® is already moderate or high 
(greater than about 0.6), the increase in droplet speed with increase of 
Moo will then play a relatively unimportant role in determining the over- 
all fractional loss-rate characteristic (on a unit distance basis). The 
reasons for this particular characteristic have been presented in some 
detail because very similar considerations apply in the cases of trends 
appearing in the remaining tabulated entries and figures • 

Additional results of the same kind are presented in table II (b) for 
an air temperature of -15° C; the remaining variables have the same 
values as in the case of the data of figure 2 (tabulated in table 11(a)). 
A detailed comparison between the two sets of results indicates that, as 
expected, the increases of fractional mass loss per unit distance which 
occur as a result of the increases of M® and M r and decreases of 1% 
are somewhat, although not substantially, more marked at the lower 
ambient temperature than at the higher. In other respects, the results 
are similar in the two cases. 

For example, at an Mr of 0 and an M r of 0.125 for the -15° C 
air -temperature case, ft changes from 0.0079 at an M^ of 0*25 to 
0.090 (a factor of 11) at an M® of 0.75; at an Mr of 0-50 and an M r 
of 0.125, ft changes from 0.00024 at an M® of 0.50 to 0-0081 (a factor 
of 34) at an M® of 0.75. 
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Figures 3 and 4 exhibit mass -loss rates as functions of Mr and 
Mjv respectively, for various sets of conditions. The data are not tab- 
ulated in the sequence in which they are plotted in these figures, but 
are included in the listings of table 11(a) . The curves of figures 2, 

3, and 4, therefore, actually present the same data in different ways. 
The various curves may be considered to be lines of intersection of 
planes with warped surfaces of three-dimensional figures. For example, 
if a surface is formed of all points for which M = 0.50, the vertical 
axis representing ft and the other axes % and M r , then the curves 
of figures 3(a) and 4(a) are lines of intersection of vertical planes 
parallel to the Mr and Mr axes, respectively. 

Interpretations of curve trends in the case of the figures 3 and 
4 are difficult; it may, however, be pointed out that, irrespective of 
time rate of loss of mass, the fractional mass-loss rates (on a unit 
distance basis) will tend toward infinity whenever both local air speed 
and relative air-droplet speed tend toward zero. For example, at Mr 
values less than 0-10, the fractional mass losses per foot of droplet 
travel exhibited in figure 3 become of the order of 0.1 as M r decreases 
below 0.125. Physically, this reflects the near absence of motion of the 
droplet. In the limit, at Mr = M r = 0, all droplets will of course 
evaporate completely. 

Droplet size effects are shown in figure 5; the same results are 
listed in table 11(c). It will be noted that while the fractional mass- 
loss rates decrease, over the droplet size range 5 to 150 microns, by a 
factor varying between 100 and 10,000, the rates of loss of actual mass 
increase, over the same droplet size range, by a factor of 270 to 2.7 
( in the same order in which the fractional loss rates are given) . 

Temperature effects are shown in figure 6- These data and a few 
additional data are tabulated in table II (d). Despite the large changes 
in vapor pressure in going from an ambient temperature of -30° C to one 
of -5° C, the factor of increase in fractional mass lost per foot varies 
between about 2 and 4 at all combinations of air-obstacle relative Mach 
number. At other altitudes, flight Mach numbers, and droplet sizes, 
this factor of increase with temperature remains substantially fixed. 

The droplet temperature increase (as is confirmed by values tab- 
ulated in this report) is usually about two thirds of the ambient tem- 
perature change; hence, the ambient vapor -pres sure increase tends to 
compensate for the liquid-surface vapor -pressure increase. 

Altitude effects are shown in figure 7 and are exhibited in numer- 
ical form in table Il(e). A remarkably small role is played by altitude 
in the determination of loss rates. 
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For example, at M w = 0.75, M r = 0, w = -25° C, and D = 15 
microns, the increase in going from 3000 to 30,000 feet is about 30 per- 
cent; at MqjS 0.75, M r = 0.25, and for the same air temperature and 
droplet diameter, the increase is about 80 percent. Under certain con- 
ditions, this change decreases and even reverses slightly. However, an 
increase of around 10 percent is typical at the higher evaporation rates. 
Mutually compensating effects are the lower relative Reynolds numbers 
(at the lower densities) on one hand and the higher mass -transfer rates 
(at the lower pressures) on the other. 

A few of the calculations (table 11(f)) were made on the basis of a 
free-stream relative humidity of 90 percent, and the balance on the basis 
of a relative humidity of 100 percent. The droplet equilibrium temper- 
atures in the latter (100 percent) case are between 0.2° and 0*5° C 
higher than in the former, while the evaporation rates are slightly 
lower, the difference usually being of the order of 2 to 7 percent. 

The results of calculations for two cases of losses from droplets 
impinging on cylinders are given in figures 8 and 9. The droplet motions 
were computed by a numerical technique; the x-motion equation alone of 
reference 3 was used, since lateral motion was assumed negligible. The 
air -flow history calculation was based on the assumption that the local 
velocity was that of incompressible potential flow, although, thereafter, 
compressible flow relations were used to determine the local air-obstacle 
relative Mach number M r , temperature, pressure, and density. When the 
x-motion equation was used, account was taken of the changing air density 
in the local relative Reynolds number. The local evaporation rate was 
obtained, for each step of the droplet motion, from the curves of fig- 
ure 3 or 4 by interpolation, the successive products of evaporation rate 
(on a unit distance basis) and step distance being then summed to obtain 
total mass losses to the respective points. 

It will be noted that the total loss in the case of the 10- 
centimeter (diameter) cylinder is less than 0.25 percent, while the 
total loss in the case of the larger (100 -centimeter) cylinder is less 
than 5 percent (the end value was obtainable only by the somewhat uncer- 
tain extrapolation indicated by the dashed area of fig. 9) . In both 
cases, most of the evaporation occurs within a distance of one cylinder 
radius from the cylinder surface. 

Finally, it is evident that dynamic calculations would yield very 
small losses. The internal time constant of a 15-micron liquid water 
sphere is 0.000023 second (fig. 10). The ratio of approach period to 
internal time constant I\j_ is 92 for the large cylinder and 9.2 for the 
smaller. It follows that, in both cases, the internal temperature remains 
substantially uniform. On the other hand, the ratio of approach period 
to external time constant is considerably smaller and varies with the 
local relative heat-transfer Nusselt number. At a distance of 2 radii 
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from the small cylinder or of about 5 radii from the large cylinder, the 
Nusselt number is not significantly greater than 2, and the ratio of 
external to internal time constant is 142, from relation (B6) . At those 
respective distances, therefore , the ratios of approach period to 
external time constant are 9.2/142 or 0*065 for the small cylinder and 
92/142 or 0-65 for the large one. At points much closer to the cylinder, 
for example, at a distance of 0.1 radius from either surface, the ratios 
are about 9.2/51 or 0-18 for the small cylinder and 92/94 for 1.0 for 
the large cylinder; the latter figures are based on relative Reynolds 
numbers of 80 and 10.8, respectively. 

It is clear that dynamic thermal response calculations based on 
external time constant alone (that is, calculations for which internal 
temperatures were considered uniform) would be valid for 15-micron 
droplet impingement on either of these cylinders. It is not necessary 
to make such calculations, however, to draw the qualitative conclusion 
that the actual total evaporation would be a few percent less than that 
given previously in the case of the larger cylinder, while the actual 
loss', in the case of the smaller one, would be considerably smaller than 
the value determined. 

It is clear that total evaporative losses from droplets approaching 
small or moderately large (up to 100 cm in the transverse direction) 
obstacles will be of the order of several percent at most for most drop- 
lets. If, however, the droplet motional history, irrespective of evap- 
oration, is such that the droplet remains for a very long (or possibly 
infinitely long) period at a finite distance from the surface, it will, 
of course, evaporate completely. Discussion of the parameters which 
determine, in particular cases, the highest droplet diameters for which 
this will occur is, however, beyond the scope of the present investiga- 
tion; such diameters, for all but the largest obstacles, will be very 
small (of the order of 5 microns). 

Duct results (quantitatively correct, to a first approximation) may 
be obtained directly from figures 2, 3, and 4 for a droplet diameter of 
15 microns, an altitude of 10,000 feet, and an ambient temperature of 
-25° C- Selection of the proper flight Mach number and stagnation Mr 
curve for the condition M r = 0 will yield the desired fractional mass 
loss rate on a unit distance basis. Selection of the M r = 0 curves is 
based on the observation that most of the duct air flow occurs at a con- 
stant ram recovery (that is, Mr value) and that, for a reasonably long 
duct and small droplet size, the droplet will be essentially at rest with 
respect to the air over most of the intra-duct droplet history. Further, 
the internal and external time constants are sufficiently small in com- 
parison with the time spent by the droplet within a duct having a length 
of at least several feet that quasi -static calculations are not seriously 
in error. 



20 


NACA TN 3024 


No droplet evaporation will occur when the duct inlet velocity- 
ratio is unity. At the other extreme, let it he assumed that the flight 
Mach number is 0.75 and the air -duct relative Mach number for intra- 
duct air is 0.25. At a droplet-air relative Mach number of 0, a value 
of £5 of 0.010 per foot is obtained from table 11(a) or figure 2(b) for 
a droplet of 15-micron diameter at 10,000-foot altitude and -25° C 
ambient temperature. Accordingly, the loss for a 10-foot duct would be 
10 percent; at a flight Mach number of unity (other conditions remaining 
the same) the loss would be about 25 percent. . Actually, the loss would 
be greater in the latter case, since account has not been taken in the 
present calculation of the increase of fractional mass loss rate with 
decrease of droplet size. It is possible to conclude that the normal 
range of droplet mass losses in ducts several feet in length is 0 to 
50 percent- 

The additional results listed in table 11(b) may also, of course, 
be used for duct or obstacle calculations. 

Factors which may be used to correct evaporation rates obtained 
from figures 2, 3, and 4 (or the data of table 11(b)) for departure of 
droplet diameter, ambient temperature, or altitude from the standard 
values (15 microns, -25° C, and 10,000 ft) assumed here may be obtained 
from figures 5, 6, and 7 or tables Il(c), 11(d), and Il(e). In each 
case, the procedure required to obtain a correction factor consists 
merely of noting the ratio of evaporation loss at the nonstandard condi- 
tion to that at the standard. For example, let the fraction of mass 
lost per foot of travel of a 10-micron droplet be required for the fol- 
lowing conditions: M^, 0.75; M^, 0.50; Mp, 0.25; ambient temperature , 

-10° C; altitude, 30,000 feet. From table 11(b) the information is 
obtained that, at 15 microns, -15° C, and 30,000 feet, the mass loss 
rate is 0.0100 per foot. From table 11(c) it can be determined that the 
ratio of loss rate at 10 microns to that at 15 microns for the given 
values of M , Mp, and M r is 19.4/10.0 = 1.94. Table Il(d) yields the 
information that the ratio of loss rate at -10° C to that at -15° C is 
12.2/10.0 = 1.22. Finally, table 11(e) yields the information that the 
ratio of loss rate at 30,000 feet to that at 10,000 feet is 7.57/6.25 = 
1.21. The product of the three factors is 2.87. The mass-loss rate 
required is therefore 0-0010X2.87 = 0.0029 per foot. The procedure is 
not rigorously valid, but the inaccuracies involved are no larger than 
the inaccuracies inherent in the basic calculations. 


CONCLUSIONS 

On the basis of a theoretical analysis of rates of loss by evapor- 
ation of atmospheric droplets moving along stagnation lines toward 
obstacles (such as wings and icing -rate measurement cylinders), it may 
be concluded that: 
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1. Little or no evaporative loss occurs from droplets approaching 
liquid-water-content measurement cylinders (that is, cylinders having 
diameters less than about 15 cm) . 

2. Evaporative losses may be as high as several percent in the case 
of small droplets approaching larger obstacles, such as wings, except 
that there is always a possibility that the droplet will never reach the 
airfoil, (in the latter case, it will, of course, evaporate completely 
if it has been approaching along the stagnation line) . 

3. Total losses from droplets moving along intake ducts, for example 

between the inlet entrance and engine screen of a jet engine, will ^ ’ 

usually be of the order of 5 to 10 percent at low temperatures for the 
smaller droplets and a fraction of that for the larger droplets, but may 
be as great as 50 percent for the smaller droplets at ambient tempera- 
tures closer to 0° C for high degrees of stagnation and moderately long 
(10 ft) ducts. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, August 17, 1953 



22 


NACA TN 3024 


APPENDIX A 
SYMBOLS 

The following symbols are used in this report: 


Ad 

droplet surface area, cur 


a o 

C3 1 (C-l - V& } q } i - *d,e,l X_1 ^ ° c2 


a l 

eg 1 (X -1 + c 2 ), °C 


Ci,c 2 ,c 3 

constants of approximate quadratic vapor -pressure equations: 


(a) (b) 



C]_ = 6136 or 5288 dyne cm ’ 2 

C 2 = 438.7 or 287.7 dyne cm ’ 2 °C’‘^ 
C 3 = 10-69 or 4.3047 dyne cm ’ 2 °C ~ 2 


c p,d 

specific heat at constant pressure of liquid water 
cal g-1 °C*1 

(droplet) , 

D 

diameter, cm 


»d 

droplet diameter, cm 


G a 

air mass -flow rate, g cm sec’ x 


h h 

heat-transfer coefficient, erg cm " 2 sec ’ 1 ^" 1 


h m 

mass-transfer coefficient, g cm ’ 2 sec ’ 1 (dyne cm’ 2 ) 

-1 

Jm 

Colburn mass -transfer parameter defined by relation 

J m = st m Sc2/3 


^f ,h 

thermal conductivity of material characteristic of 
erg cm -1 °C _1 sec -1 

film, 

kf ,ra 

mass-transfer conductivity of material characteristic of film 

P a 3 /Md\ 0 n 

and of diffusing material, — — l jjj-), ® cra - ^ sec" x 

(dyne cm“^) -G 

X h,a 

thermal conductivity of air, erg cm ’ 1 °C ’ 1 sec ’ 1 
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k h,d 

% 

Md 

Mr 

Mr 


N.P. 

Wu h 

NUjjj 

P 

Pa 

P e,l,r 

P st,® 

Pr 

^d, s 

p d,co 

P d,e, Z 

^m 

R 

Re 


thermal conductivity of droplet, erg cm“P °C”P sec"^ 

3 3-1 

molecular weight of dry air, g (22.4X10 cm )” 
molecular weight of water, g (22.4X10^ cm^)“P 
local Mach number of air flow with respect to obstacle 
local relative air -droplet Mach number 
flight Mach number 

(1) constant equal to 1, 2, or 3 as discussed in text in 
appendix B 

(2) mass -transfer rate, g sec"P cm -2 
normal pressure (l atm), 1,013,250 dyne cm -2 
Nusselt number for heat transfer, h^D/k^ ^ 

Nusselt number for mass transfer, h^^/k^ m 
air pressure, dyne cm -2 

o 

mean local dry-air partial pressure, dyne cm 
local effective relative air pressure, dyne cm -2 

static ambient flight air pressure (at infinity), dyne cm -2 

Prandtl number 

droplet equilibrium (psychrometric) vapor pressure, dyne cm - 

ambient static flight water vapor pressure (at infinity), 
dyne cm* 2 

local effective water vapor pressure, dyne cm -2 

rate of loss of mass of droplet, g sec"P 
principal (ref.) dimension of obstacle, cm 
Reynolds number 
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d,e,Z 


d,t,l 


L St, 09 


"d,e,Z 


St, CD 


relative air-droplet Reynolds number based upon local static 
air density, droplet diameter, relative air -droplet air 
speed, and viscosity at the relative air-droplet total 
temperature 

droplet radius, cm 

Schmidt number, \i / p (3 

SI cl 

Stanton number for mass tr ans fer 
absolute temperature, °K 
volume-mean droplet temperature, °K 

effective local air temperature with respect to droplet, °K 

air total temperature with respect to droplet, °K 

local static air temperature, °K 

ambient (flight) static air temperature, °K 

air total tenperature with respect to obstacle, °K 
temperature, °C 

effective local air temperature with respect to droplet, °C 
equilibrium (psychrometric) droplet temperature, °C 
ambient (flight) static air temperature, °C 

local air speed with respect to obstacle, cm sec -1 

f light obstacle speed, cm sec -2- 

local air -droplet relative speed, cm sec -1 

* h n/ h h> ° c ( d yne cm' 2 )' 1 

diffusion coefficient, cm 2 " sec -2- 

diffusion coefficient at 0® C and atmospheric pressure, 
cm* 2 sec -2 - 
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r ^ number of internal sphere time constants per unit of time 

required for obstacle to move through air at speed a 

distance R 

X ratio of specific heats 

6 time, sec 

ext external time constant of object (sphere), sec 

®d,int internal time constant of object (sphere), sec 

thermal diffusivity of droplet (water), cm** sec - ^- 
A time rate of loss of fractional mass, sec - -*- 

X heat of vaporization, erg g - -*- 

q viscosity of air, poise 

viscosity of air at temperature t, l > poise 

p density, g cm“3 

droplet (water) density, g cm - ^ 
p local static air density, g cm - ^ 

St y L 


a flight relative humidity 

0 fractional droplet mass lost per centimeter of droplet travel 

toward obstacle; loss per foot is given in tabulations and 
graphs 

Subscripts : 

a appertaining to air 

d appertaining to droplet 

e effective value 

f appertaining to film 

h heat transfer 

1 local 
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m mass transfer 

N-P* normal (atmospheric) pressure 

N.T.P. normal (atmospheric) pressure and (0° C) temperature 
R air-obstacle relative value 

r droplet -air relative value 

r- 

s equilibrium (psychrometric) value 0 

st static value 

t total value 

« value in undisturbed stream 
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APPENDIX B 


SPHERE TIME CONSTANTS 

The thermal response of a solid body may, to a first approximation, 
be described in terms of what may be designated an "internal time con- 
stant" and an "external time constant." (Although liquid spheres are 
discussed in this report, they are treated as though solid.) For pre- 
sent purposes the internal time constant is designated 6 ^ int and is 

defined as the time required for the mean sphere temperature to change 
by l-e~l of the asymptotic change upon imposition of an instantaneous 
surface temperature change. The external time constant is designated 
%,ext and is defined as the time required for the sphere temperature 

(assumed uniform throughout) to change by 1-e" 1 of the asymptotic 
change upon imposition of an instantaneous environmental temperature 
change. 

From curve V, figure 9, page 84 of reference 15, it is found that 

O 

the Fourier modulus K^O/r^ has, in the case of the sphere, the value 

0.056 when the temperature ratio T^/T s = 0.632; consequently, the 

internal time constant of a sphere is given approximately by the 
relation 


%, int “ °' 056 r d/ K d ( B1 ) 

In the case of water, K. = 0.00136 at 0° C, in c.g.s. -calorie units, 

^ 2 7 

and it is found that e d ^ lnt = 41. r d ; r d is here given in centimeters. 
Figure 10 exhibits this relation over the range D d = 2r d = 1 to 1000 
microns, over which 0 d ,i n t varies between ~10 -7 and 10" 1 second. At 
a typical droplet size, 10 microns, the time constant is about 10 -5 
second. 

The question arises as to the magnitude of this period as compared 
with the approach time of a droplet. 

If i S defined as the number of internal sphere time constants 

per unit of time required for the air (or droplet) at infinity to move 
(relatively) through a distance equal to the selected principal dimen- 
sion of the obstacle, the relation 
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applies (the cylinder radius has "been selected as the principal dimen- 
sion, rather than the diameter). 

Values of T i are plotted in figure 11 as a function of U w for 
a number of combinations of R and D d . The range of r d covered is 
about 0.01 to 1000, of U<» 225 to 825 miles per hour, of R 0-5 to 
180 centimeters, and of D d 6-25 to 400 microns. It is clear that for 
large (>5, say) values of the sphere temperature will be uniform 

except possibly in regions of maximum rate of air temperature change, 
where small gradients will exist. On the other hand, it may be stated 
without proof that for small (<0.2) values, it is possible to make an 
estimate of the fractional volume (adjacent to the surface) of the 
sphere which alone participates, to a first approximation, in the ther- 
mal exchanges with the environment, and so to make rough calculations 
using, in effect, a hollow sphere in the liquid portion of which the 
temperature is taken to be uniform. In the region 0-2 < 1^ <5, roughly, 
no simple approach appears to be possible; this region is accordingly 
labeled "critical region" in figure 2. 

It will be noted that many combinations of liquid-water -content 
measuring cylinder diameter and droplet diameter fall within or near the 
critical region. While it is possible to deduce that in such cases the 
droplet surface temperatures will be higher than they will be for first- 
order response (that is, the response of the sphere as a unit), it is not 
easy to make quantitative evaluations. 

In the subcritical region the (approximate) calculations based on 
the "shell" concept, while possibly easier, will lead to results varying- 
very widely with the particular conditions, and such calculations are 
accordingly not made in the present report. Again, in general, surface 
temperatures will be higher, and usually markedly so, than those attained 
on. the basis of first-order response. 


The external time constant %^ ex t in the case of heat transfer 

without evaporation appears, for example, as a by-product of investiga- 
tions of heat transfer to and from solid bodies through uniform boundary 
layers as carried out by several investigators. Jakob (ref. 16, pp. 
270-291) summarizes their results and, in particular, presents (p. 291) 
a chart in which exponential changes of (uniform) plate, cylinder and 

mfa h h 0 

sphere temperatures are given as functions of the variable — — — — . in 


k h,d r d 

that expression, m has the values 1, 2, and 3 for plate, cylinder, and 
sphere, respectively, and r d is half the plate thickness or, otherwise, 
the body radius. When the variable in question equals unity, 0 equals 
®d, ext b y definition and, in the case of the sphere, the relation 
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0 - v 

d, ext " 3kd hh 


d r d 


p d c p,d r d 
3hv, 


(B3) 


is obtained. The same relation is, of course, obtainable very directly 
and simply by considering that the system sphere— boundary layer behaves 
like anelectrical resistance -capacity circuit of capacity 

(4*/3)r d p d c p d and resistance - — V— , the product RC then yield- 

4n r£ h h 

ing % ext. For “»ny Purposes it is preferable to use the Russelt num- 
ber of heat transfer 


Nuu = 


h h D d 

^h,a 


when Nu^ is used, this equation becomes 


0 


d, ext “ 3k 


2r, p c , 
d d p,d 


h,a 


Nu v 


(B4) 


It is known that 


Nu h 


0.6 

OC Re r 


and it is instructive to note, therefore, that 


®d, ext 


OC 



Pd c p,d 

kh,a 


(Obviously, air relative speed, air density, and air viscosity factors 
have teen omitted. ) 


In view of the fact that, in contra-distinction to 0, . . , Q, 

d,int> d, ex 

depends on a comparatively large number of factors, it did not appear 
worthwhile to present any set of curves from which its magnitude might 
be estimated. On the other hand, a direct comparison with 0-, is 

fruitful. It is clear that d, int 


0 


d,ext 


11.9 k. 


h,d 


®d, int kh, a Nu^ 


(B5) 


Upon evaluation of the thermal conductivities of air and liquid water ' 
at 0° C, it is found that 


^d, ext _ 284 
%, int = 


(B6) 


in the case of a sphere of liquid water in air at 0° C* 
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It is generally accepted, on the basis of both theory and observa- 
tion (ref. 11, for example), that the minimum value of N% is 2. The 
maximum value that the relative Reynolds number of a droplet-air system 
is likely to reach, for flight Mach numbers less than unity, is about 
1000, at which Reynolds number the value of Nuh (ref. 11, fig. 9) will 
be about 20. The range of the ratio of equation (B6) will therefore be 
about 10:1, from values of about 140 (extreme high) through typical 
values of the order of 40 to values of about 14 (extreme low). 

When evaporation (or, in general, mass transfer) is occurring, the 
situation is not so simple since latent heat of vaporization or conden- 
sation may be either hastening or slowing down a particular droplet 
temperature change. It is to be noted that if the environment of a 
large droplet is suddenly altered by an adiabatic compression, some con- 
densation may conceivably occur before the instant at which the droplet- 
surface vapor pressure, as a result of surface -temperature rise, once 
again equals or exceeds the increased ambient water vapor pressure. 
Initially, then, the addition of heat of condensation to droplet heat 
content would tend to accelerate surface temperature rise and, therefore, 
decrease the initial time constant. 

In general, then, even the first-order response cannot be calcu- 
lated unless assumptions are made concerning condensation rates in the 
case of larger droplets or unless the situation is that of approach of 
small droplets (5 to 20 microns) to obstacles that are not small. 
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TABLE I. - COMPARISON OF ACTUAL VAPOR PRESSURES AND PRESSURES 
GIVEN BY EQUATIONS (2a) AND (2b) 



Temper- 

ature, 

°C 

Actual 
vapor 
pressure, 
dyne cm -2 

Pressure 
(eq. (2a)) 

Error , 
percent 
(eq. (2a)) 

Pressure 
(eq. (2b)) 

Error , 
percent 
(eq-. (2b)) 

-32.0 

421 



490 

+16.4 

-31.5 

441 



497 

+12.7 

-31.0 

463 



506 

+9.3 

-30.5 

485 



518 

+6.8 

-30.0 

509 



531 

+4.3 

-29.5 

533 



547 

+2.6 

-29 .0 

559 



565 

+1.1 

-28.5 

586 



585 

-.17 

-28.0 

613 



607 

-.98 

-27.5 

643 



632 

-1.7 

1 

ro 

o 

673 



658 

-2.2 

-26.5 

704 



687 

-2.4 

-26.0 

737 



718 

-2.6 

-25.5 

771 



751 

-2.6 

-25.0 

807 



786 

-2.6 

-24.5 

844 



823 

-2.5 

-24.0 

883 



863 

-2.3 

-23.5 

923 



905 

-2.0 

-23.0 

965 



948 

-1.8 

-22.5 

1008 



994 

-1.4 

-22.0 

1054 



1042 

-1.1 

-21.5 

1101 



1093 

-.73 

-21.0 

1150 



1144 

-.52 

-20.5 

1201 



1199 

-.17 

-20.0 

1254 



1256 

+ .16 

-19.5 

1309 



1315 

+ .46 

-19.0 

1366 



1376 

+ .73 

-18.5 

1426 



1439 

+ .91 

-18.0 

1488 



1504 

+1.1 

-17.5 

1552 

1733 

+11.7 

1572 

+1.3 

-17.0 

1619 

1768 

+9.2 

1641 

+1.4 
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TABLE I. - Continued. COMPARISON OF ACTUAL VAPOR PRESSURES AND 
PRESSURES GIVEN BY EQUATIONS (2a) AND (2b) 


Temper- 

ature, 

°C 

Actual 
vapor 
pressure, 
dyne cm” 

Pressure 
(eq. (2a)) 

Error, 
percent 
(eq. (2a)) 

Pressure 
(eq. (2b)) 

Error , 
percent 
(eq. (2b)) 

-16.5 

1688 

1808 

+7.1 

1713 

+1.5 

-16.0 

1760 

1853 

+5.3 

1787’ 

+1.5 

-15.5 

1834 

1905 

+3.9 

1863 

+1.6 

-15.0 

1912 

1961 

+2.6 

1941 

+1.5 

-14.5 

1992 

2023 

+1.6 

2022 

+1.5 

-14.0 

2076 

2089 

+ .63 

2104 

+1.3 

-13.5 

2162 

2162 

0 

2189 

+1.2 

-13.0 

2252 

2239 

-.50 

2275 

+1.0 

-12.5 

2345 

2323 

-.94 

2365 

+ .85 

-12.0 

2441 

2411 

-1.2 

2456 

+ .61 

-11.5 

2541 

2505 

-1.4 

2549 

+ .31 

-11.0 

2644 

2604 

-1.5 

2644 

0 

-10.5 

2752 

2709 

-1.6 

2742 

-.36 

-10.0 

2863 

2818 

-1.6 

2841 

-.77 

-9.5 

2978 

2933 

--1.5 

2943 

-1.2 

-9.0 

3097 

3054 

-1.4 

3048 

-1.6 

-8.5 

3221 

3179 

-1.3 

3154 

-2.1 

-8.0 

3348 

3311 

-1.1 

3262 

-2.6 

- 7.5 

3481 

3447 

-.98 

3372 

-3.1 

-7.0 

3618 

3589 

-.80 

3485 

-3.7 

-6.5 

3759 

3736 

-.61 

3600 

-4.2 

-6.0 

3906 

3889 

-.44 

3717 

-4.8 

-5.5 

4058 

4047 

-.27 



-5.0 

4215 

4210 

-.12 



-4.5 

4377 

4378 

+ .12 



-4.0 

4545 

4552 

+ .15 



-3.5 

4719 

4732 

+ .28 



-3.0 

4898 

4916 

+ .37 



-2.5 

5084 

5106 | 

+ .43 



-2.0 

5275 

5301 

+ .49 



-1.5 

5473 

5502 

+ .53 



-1.0 

5678 

5708 

+ .53 
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TABLE I.- Concluded. COMPARISON OF ACTUAL VAPOR PRESSURES AND 
PRESSURES GIVEN BY EQUATIONS (2a) AND (2b) 



Temper- 

ature, 

OC 

Actual 
vapor 
pressure, 
dyne cm" 2 

Pressure 
(eg. (2a)) 

Error, 
percent 
(eg. (2a)) 

Pressure 
(eg. (2b)) 

Error, 
percent 
(eg. (2b)) 

-0.5 

.5889 

5919 

+ .50 



0 

6108 

6136 

+ .46 



.5 

6333 

6358 

+ .39 



1.0 

6566 

6585 

+ .29 



1.5 

6807 

6818 

+ .16 



2.0 

7055 

7056 

+ .01 



2.5 

7311 

7300 

-.15 



3.0 

7575 

7548 

-.36 



3.5 

7848 

7802 

-.59 



4.0 

8129 

8062 

-.82 



4.5 

8420 

8327 

-1.1 



5.0 

8719 

8597 

-1.4 



5.5 

9028 

8872 

-1.7 



6.0 

9347 

9153 

-2.1 



6.5 

9675 

9439 

-2.4 



7.0 

10013 

9731 

-2.8 



7.5 

10362 

10028 

-3.2. 



8.0 

10722 

10330 

-3.7 



8.5 

11092 

10637 

-4.1 



9.0 

11474 

10950 

-4.6 



9.5 

11867 

11268 

-5 0 



10.0 

12272 

11592 

-5.5 



10.5 

12690 

11921 

-6.1 



11.0 

13119 

12255 

-6.6 



11.5 

13562 

12595 

-7.1 



12.0 

14017 

12940 

-7.7 



12.5 

14486 

13291 

-8.2 



13.0 

14969 

13646 

-8.8 



13.5 

15466 

14007 

-9.4 



14.0 

15977 

14373 

-10.0 



14.5 

16503 

14745 

-10.7 



15.0 

17044 

15122 

-11.3 
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TABLE II - RESULTS 


(a) Flight Mach number effects at -25 


Altitude, t 
ft) 


V r * 

t/sec 

Re r 

^,8’ 

°C 

65.68 

18.12 

-22.29 

66.86 

19.64 

-14.94 

68.83 

22.32 

-3.86 

71.52 

26.51 

9.66 

130.3 

35.88 

-21.88 

132.7 

38.90 

-14.52 

136.6 

44.28 

-3.44 

141.9 

52.46 

10.14 

260.6 

71.30 

-20.33 

265.4 

77.25 

-12.96 

273.1 

87.95 

-1.96 

283.8 

104.3 

11.79 

530.8 

150.3 

-6.98 

546.3 

171.0 

3.69 

567.6 

202.8 

a 17 .89 

819.5 

245.9 

12.87 

851.4 

291.8 

a 28 .01 

0 

0 

-17.45 

-5.98 

7.78 

65.26 

17.64 

-24.88 

66.50 

19.14 

-17.42 

68.40 

21.74 

-6.11 

71.09 

25.79 

7.44 

129.5 

34.92 

-24.46 

131.9 

37.83 

-17.03 

135.8 

43.05 

-5.71 

141.1 

51.03 

7.90 

259.0 

69.38 

-22.89 

263.9 

75.17 

-15.47 

271.5 

85.50 

-4.21 

282.2 

101.4 

9.52 

527.9 

146.2 

-9.38 

543.0 

166.4 

1.42 

564.3 

197.3 

a 15.50 

814.7 

239.0 

10.52 

846.5 

283.7 

a 25.43 

0 

0 

-13.11 

1.19 

65.26 

17.64 

-24.88 

67.16 

20.03 

-13.13 

69.78 

23.75 

.96 

129.5 

34.92 

-24.46 

133.2 

39.69 

-12.73 

138.5 

47.01 

1.38 

259.0 

69.38 

-22.89 

266.5 

78.83 

-11.14 

276.9 

93.44 

2.91 

518.1 

134.8 

-16.77 

533.2 

153.3 

-5.29 

553.8 

181.7 

8.66 

799.6 

220.2 

3.70 

830.7 

261.2 

a 18 .09 

0 

0 

-9.49 

65.26 

17.64 

-24.88 

67.78 

20.85 

-9.57 

129.5 

34.92 

-24.47 

134.5 

41 .31 

-9.14 

259.0 

69.38 

-22.89 

269.0 

82.10 

-7.63 

518.1 

134.8 

-16.77 

538.1 

159.5 

-1.97 

776.9 

193.4 

-7.12 

807.1 

229.2 

7.03 


1.41x10 

4.98 


A, Ui + V r , 
fractional ft/sec 
mass/sec 


23.6 1333 

51.0 1384 


‘Droplet temperatures emu evapwauwi 1 « 

Evaporation rates omitted for droplet temperatures greater than 25° C. 
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TABLE II - Continued. RESULTS 
(c) Droplet diameter effects, 


Altitude, 
ft °i 


fractional ft/sec 
mass/sec 


10,000 -15 0.75 0 


4.57X10' 

10.6 

17.6 

34.2 

87.3 
418 


.75 835.6 79.20 a f20.01 

158.4 20.35 

237.6 20.58 

396.0 20.92 

792.0 21.39 

2376 22.22 


1.94X10' 

5.20 

9.38 

20.1 

57.3 
313 

2.66x10’ 

5.32 

7.99 

13.3 
26.6 
79.9 


1.75X10 

4.66 

8.43 

18.0 

51.3 

276 

1 .39x10”' 
2.79 
4.18 
6.97 
13.9 
41.8 


.332x10’ 

.828 

1.44 

2.96 

8.06 

41.8 


a Droplet temperatures and evaporation rates not quantitative for calculated droplet temperatures higher than about 15° C. 
Evaporation rates omitted for droplet temperatures greater than 25° C. 
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TABLE II - Continued. RESULTS 


(d) Temperature effects. 


Dd- 

microns 

Altitude, 

ft 

°C 

M <D 

Mr 

M r 

v r * 

ft/sec 

Re r 

, s’ 
°c 

Qm' 

g/sec 

A, 

fractional 

mass/sec 

Ui + v r , 
ft/sec 

n, 

fractional 

mass/ft 

15 

10,000 

-30 

0.75 

0 

0.063 

68.14 

22.76 

-7.75 

8.63X10 -9 

4.88 

68.14 

71.6X10 -3 



-25 




68.83 

22.32 

-3.86 

11.5 

6.53 

68.83 

94.9 



-20 




69.52 

21.93 

.03 

14.5 

8.22 

69.52 

118 



-15 




70.18 

21.57 

3.93 

17.6 

9.98 

70.18 

142 



-10 




70.90 

21.23 

7.82 

20.8 

11.8 

70.90 

166 



-5 




71.59 

20.93 

12.04 

23.5 

13.3 

71.59 

186 



-30 



.25 

270.3 

89.58 

-5.78 

15.8 

8.94 

270.3 

33.1 



-25 




273.1 

87.95 

-1.96 

20.8 

11.8 

273.1 

43.2 



-20 




275.9 

. 86.36 

1.92 

25.9 

14.7 

275.9 

53.2 



-15 




278.5 

84.95 

5.79 

31.1 

17.6 

278.5 

63.2 



-10 




281.3 

83.63 

9.69 

36.4 

20.6 

281.3 

73.2 



-5 




284.0 

82.44 

13.91 

40.8 

23.1 

284.0 

81.4 



-30 



.75 

811.0 

250.2 

9.18 

5.89X10~ 8 

33.3 

811.0 

41.1 



-25 




819.5 

245.9 

„ 12.87 

7.09 

40.1 

019.5 

48.9 



-20 




827.7 

241.7 


16.68 

8.25 

46.7 

827.7 

56.4 



-15 




835.6 

237.6 


20.58 

9.38 

53.1 

835.6 

63.5 



-10 




844.1 

233.9 


24.59 

10.5 

59.5 

844.1 

70.5 



-5 




852.0 

230.7 


120.96 



852.0 




-30 


.25 

0 

0 

0 


-10.09 

3.66X10 -9 

2.07 

268.8 

7.70 



-25 






-5.98 

4.98 

2.82 

271.5 

10.4 



-20 






-1.96 

6.43 

3.64 

274.3 

13.3 



-15 






2.05 

7.99 

4.52 

276.9 

16.3 



-10 






6.01 

9.60 

5.43 

279.6 

19.4 



-5 






10.26 

11.0 

6.22 

282.2 

22.1 



-30 



.063 

67.75 

22.16 

-10.18 

7.00 

3.96 

336.6 

11.8 



-25 




68.40 

21.74 

-6.11 

9.43 

5.34 

339.9 

15.7 



-20 




69.13 

21.36 

-2.12 

12.1 

6.83 

343.4 

19.9 



-15 




69.78 

20.99 

1.84 

14.0 

8.39 

346.7 

24.2 



-10 




70.47 

20.66 

5.77 

17.7 

10.0 

350.2 

20.6 



-5 




71.13 

20.36 

10.02 

20.1 

11.4 

353.4 

32.3 



-30 



.25 

268.8 

87.20 

-8.11 

12.6 

7.15 

537.7 

13.3 



-25 




271.5 

85.50 ! 

-4.21 

17.2 

9.73 

543.0 

17.9 



-20 




274.3 

84.00 

-.27 

21.7 

12.3 

548.5 

22.4 



-15 




276.9 

82.59 

3.66 

26.4 

14.9 

553.8 

27.0 



-10 




279.6 

81.26 

7.59 

31.2 

17.7 

559.4 

31.6 



-5 




282.2 

80.15 

11.83 

35.3 

20.0 

564.6 

35.4 



-30 



.75 

806.4 

243.5 

6.80 

51.6 

29.2 

1075 

27.1 



-25 




814.7 

239.0 

10.52 

62.7 

35.5 

1086 

32.7 



-20 




822.8 

234.8 


14.33 

73.6 

41.6 

1097 

38.0 



-15 




830.7 

230.9 


18.23 

84.3 

47.7 

1108 

43.0 



-10 




838.9 

227.3 


22.21 

94.8 

53.6 

1110 

48.0 



-5 




846.8 

224.1 


_26 . 58 



1129 




-30 


.50 

0 

0 

0 

-17.56 

1.76 

.995 

527.9 

1.08 



-25 






-13.11 

2.52 

1.43 

533.1 

2.68 



-20 






-8.50 

3.19 

1.81 

538.4 

3.35 



-15 






-4.24 

4.18 

2.36 

543.6 

4 .35 



-10 






-.07 

5.31 

3.00 

549.2 

5.47 



-5 






4.32 

6.28 

3.55 

554.1 

6.41 



-30 



.063 

66.50 

20.43 


-17.56 

3.31 

1.87 

594.4 

3.15 



-25 




67.16 

20.03 


-13.13 

4.71 

2.67 

600.3 

4.44 



-20 




67.85 

19.68 

-8.55 

5.91 

3.34 

606.6 

5.51 



-15 




68.50 

19.32 

-4.31 

7.66 

4.34 

612.1 

7.08 



-10 




69.19 

19.04 

-.16 

9.65 

5.46 

618.4 

0.83 



-5 




69.85 

18.77 

4.20 

11.3 

6.39 

624.3 

10.2 



-30 



.25 

263.9 

80.39 

-15.54 

6.41 

3.63 

791.7 

4.58 



-25 




266.5 

78.83 


-11.14 

8.83 

5.00 

799.5 

6.25 



-20 




269.2 

77.48 

-6.75 

11.3 

6.41 

807.7 

7.94 



-15 




271.8 

76.10 

-2.55 

14.4 

8.17 

015.6 

10.0 



-10 




274.6 

74.90 

1.53 

17.0 

10.1 

823.8 

12.2 



-5 




277.1 

73.85 

5.87 

20.7 

11.7 

831.4 

14.1 



-30 



.75 

791.7 

224.3 

-.08 

32.9 

18.9 

1320 

14.1 



-25 




799.6 

220.2 

3.70 

41.8 

23.6 

1333 

17.7 



-20 




807.7 

216.3 

7.56 

50-. 5 

20.6 

1346 

21.2 



-15 




815.6 

212.7 

„ 11.48 

59.2 

33.5 

1359 

24.6 



-10 




823.8 

209.4 

* 

15.45 

67.8 

38.4 

1373 

27.9 



-5 




831.4 

206.6 


JL9.77 

75.1 

42.5 

1385 

30.7 


a Droplet temperatures and evaporation rates not quantitative for calculated droplet temperatures higher than about 15° C. 
Evaporation rates omitted for droplet temperatures greater than 25° C. 
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TABLE II - Continued. RESULTS 


(e) Altitude effects. 


microns 

Altitude, 

ft 

°c 


«R 

«r 

V 

ft/sec 

Re r 

^,3- 

°C 

g/sec 

A, 

fractional 

mass/sec 

U l + v r , 
ft/sec 

0, 

fractional 

mass/ft 

15 

3,000 

-25 

0.75 

0 

0.063 

68.83 

29.58 

-2.74 

10.6X10* 9 

5.98 

68.83 

86.8X10' 3 


6,000 






26.21 

-3-22 

11.0 

6.21 


90.3 


10,000 






22.32 

-3.86 

11.5 

6.53 


94.9 


15,000 






10.53 

-4.64 

12.2 

6.91 


100 


22, S00 






14.11 

-5.83 

13.2 

7.44 


100 


30,000 






11.08 

-6.93 

14.0 

7.90 


115 


3,000 




.25 

273.1 

116.4 

-.70 

19.7 

11.2 

273.1 

40.8 


6,000 






103.2 

-1.24 

20.2 

11.4 


41.8 


10,000 






87.93 

-1.96 

20.8 

11.0 


43.2 


15,000 






72.95 

-2.83 

21.5 

12.2 


44.6 


22,500 






55.56 

-4.17 

22.4 

12.7 


46.4 


30,000 






43.64 

-5.39 

23.0 

13.0 


47.7 


3,000 




.75 

819.5 

325.7 

a l5. 16 

69.3 

39.2 


47.8 


6,000 






288.3 

14.19 

70.0 

39.6 


48.4 


10,000 






245.9 

12.87 

70.9 

40.1 


48.9 


15,000 






204.0 

11.29 

71.5 

40.5 


49.4 


22,500 






155.3 

8.910 

71.6 

40.5 


49.5 


30,000 






122.0 

6.805 

71.3 

40.3 


49.2 


3,000 



.25 

0 

0 

0 

-5.06 

4.16 

2.35 

271.5 

8.66 


6,000 







-5.45 

4.50 

2.55 


9.38 


10,000 







-5.98 

4.98 

2.82 


10.4 


15,000 







-6.64 

5.58 

3.15 


11.6 


22,500 







-7.65 

6.51 

3.68 


13.6 


30,000 







-8.56 

7.33 

4.15 


15.3 


3,000 




.063 

68.40 

28.80 

-5.15 

8.55 

4.84 

339.9 

14.2 


6,000 






25.52 

-5.56 

0.92 

5.05 


14.9 


10,000 






21.74 

-6.11 

9.43 

5.34 


15.7 


15,000 






18.03 

-6.60 

10.0 

5.66 


16.7 


22,500 






13 . 74 

-7.85 

10.9 

6.15 


18.1 


30,000 






10.79 

-0.79 

11.6 

6.56 


19.3 


3,000 




.25 

271.5 

113.3 

-3.14 

16.1 

9.12 

543.0 

16.8 


6,000 






100.4 

-3.57 

16.6 

9.40 


17.3 


10,000 






85.50 

-4.21 

17.2 

9.73 


17.9 


15,000 






70.95 

-4.99 

17.8 

10.1 


18.6 


22,500 






54.05 

-6.16 

18.7 

10.6 


19.5 


30,000 






42.45 

-7.26 

19.3 

10.9 


20.1 


3,000 




.75 

814.7 

316.7 

12.61 

61.1 

34.6 

1086 

31.8 


6,000 






280.7 

11.73 

61.9 

35.0 


32.3 


10,000 






239.0 

10.52 

62.7 

35.5 


32.7 


15,000 






190.3 

9.04 

63.3 

35.8 


33.0 


22, 500 






151.1 

6.86 

63.7 

36.0 


33.2 


30,000 






118.7 

4.90 

63.5 

35.9 


33.1 


3,000 



.50 

0 

0 

0 j 

-12.57 

2.05 

1.16 

533.1 

2.17 


6,000 







-12.00 

2.24 

1.27 


2.38 


10,000 







-13.11 

2.52 

1.43 


2.68 


15,000 







-13.51 

2.80 

1.63 


3.05 


22,500 







-14.17 

3.43 

1.94 


3.64 


30,000 







-14.78 

3.98 

2.25 


4.23 


3,000 




.063 

67.16 

26.52 

-12.63 

4.11 

2.32 

600.3 

3.87 


6,000 






23.51 

-12.80 

4.38 

2.48 


4.13 


10,000 






20.03 

-13.13 

4.71 

2.67 


4.44 


15,000 






16.62 

-13.57 

5.10 

2.89 


4.81 


22,500 






12.66 

-14 . 23 

5.70 

3.23 


5.38 


30,000 






9.939 

-14.89 

6.21 

3.51 


5.85 


3,000 




.25 

266.5 

104.4 

-10.46 

8.02 

4.54 

799.5 

5.67 


6,000 






92.43 

-10.75 

8.37 

4.73 


5.92 


10,000 






78.63 

-11.14 

8.83 

5.00 


6.25 


15,000 






65.37 

-11.65 

9.35 

5.29 


6.62 


22,500 






49.80 

-12.45 

10.1 

5.72 


7.15 


30,000 






39.12 

-13.23 

10.7 

6.05 


7.57 


3,000 




.75 

799.6 

291.5 

5.29 

40.3 

22.8 

1333 

17.1 


6,000 






250.2 

4.61 

41.0 

23.2 


17.4 


10,000 






220.2 

3.70 

41.8 

23.6 


17.7 


15,000 






182.6 

2.58 

42.4 

24.0 


18.0 


22,500 






139.1 

.91 

43.1 

24.4 


10.3 


30,000 






109.2 

-.61 

43.4 

24.5 


18.4 


a Broplet temperatures and evaporation rates not quantitative for calculated droplet temperatures higher than about 15° C. 
Evaporation rates omitted for droplet temperatures greater than 25° C. 
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TABLE III. - AMBIENT CONDITIONS OF CALCULATIONS 

[Pressures from Signal Corp. RAOB chart (ML-248) 
for particular temperatures and altitudes^ 


Altitude, 

Flight air temperature, 

* Oz-t 

Air pressure 

ft 


dyne cm"^ 

atm 

3,000 

-30 

8.925X10 5 

0.8808 


-25 

8.94 

.8823 


-20 

8.96 

.8843 


-15 

8.98 

.8863 


-10 

8.995 

.8877 


-5 

9.015 

.8897 

6,000 

-30 

7.89 

0.7787 


-25 

7.92 

.7816 


-20 

7.955 

.7851 


-15 

7.995 

.7890 


-10 

8.025 

.7920 


-5 

8.06 

.7955 

10,000 

-30 

6.70 

0.6612 


-25 

6.75 

.6662 


-20 

6.80 

• 6711 


-15 

6.85 

. 6760 


-10 

6.905 

.6815 


-5 

6.97 

.6879- 

15,000 

-30 

5.54 

0.5468 


-25 

5.60 

.5528 


-20 

5.66 

.5586 


-15 

5.72 

.5645 


-10 

5.79 

.5714 


-5 

5.86 

.5783 

22,500 

-30 

4.185 

0.4130 


-25 

4.265 

.4209 


-20 

4.345 

.4288 


-15 

4.42 

.4362 


-10 

4.50 

.4441 


-5 

4.59 

.4530 

30,000 

-30 

3.26 

0.3217 


-25 

3.35 

.3306 


-20 

3.44 

.3395 


-15 

3.53 

.3484 


-10 

3.625 

.3578 


-5 

3.725 

.3676 
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Figure 1. - Motional relatione among air-stream, droplet, and obstacle. 
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.001 1 1 1 1 1 1 1 1 1 1 1 i 1 i 1 i i 

(a) Local relative air-obstacle Mach number, Mr, 0. 



Flight Mach number, M^ 

(b) Local relative air-obstacle Mach number, (d) Local relative air-obstacle Mach 

Mr, 0.25. ' number, Mr, 0.75. 

Figure 2. - Variation of mass loss rate with flight Mach number. Altitude, 10,000 feet; ambient tenders ture , 
-25° C; droplet diameter, 15 microns. 
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l 

! 


0 .1 .2 .3 .4 .5 .6 


Local relative air-obstacle Mach number, Mr 


(a) Flight Mach number, M^, 0.50. 


Figure 3. - Variation of mass loss rate with relative air-obstacle Mach number 
Altitude, 10,000 feet; ambient temperature, -25° C; droplet diameter, 

15 microns . 
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Local relative air- obstacle Mach number, 

(c) Flight Mach number, M^, 1.00. 

Figure 3. - Concluded. Variation of mass loss rate with 
relative air-obstacle Mach number. Altitude, 10,000 feet; 
ambient temperature, -25° C; droplet diameter, 15 microns. 
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NACA TN 3024 



Figure 4. - Variation of mass loss rate with relative air-droplet Mach number. 
Altitude, 10,000 feet; ambient temperature, 25° C; droplet diameter, 

15 microns. 
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Figure 5. - Variation of mass-loss rate with droplet diameter. Altitude, 10,000 feet; ambient 
temperature, 15° C; flight Mach number, 0.75. 
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(b) Local relative air-obstacle Mach number, Mr, 0.25. 

Figure 5. - Continued. Variation of mass-loss rate with droplet diameter. Altitude, 10,000 feet; 
ambient temperature, 15° C; flight Mach number, 0.75. 
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-obstacle relative 



Figure 8. - Air and droplet speeds along stagnation lines ahead of circular cylinders. Altitude, 
10,000 feet; ambient temperature, -25° C; droplet diameter, 15 microns; flight Mach number, 0.75 
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NACA TN 5024 



Figure 9. - Limiting mass losses from droplets impinging on circular cylinders. 
Altitude, 10,000 feet; ambient temperature, -25° C; droplet diameter, 15 microns; 
flight Mach number, 0.75. 









